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Abstract 

CM . 

£NJ , We investigate the form factors for exclusive semileptonic decays of i?-meson to D, D* 

based on the parton picture and helped by the results of the HQET. We obtain the numerical 
results for the slope of the Isgur-Wise function, which is consistent with the experimental re- 
sults, and we extracte the dependences of hadronic form factors on q 2 by varying input heavy 
qq \ quark fragmentation function without the nearest pole dominance ansatze. 
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*th . 1 Introduction 

a: 

Z?-meson decay processes have been studied in detail as providing many interesting informations 
on the interplay of electroweak and strong interactions and as a source extracting the parameters of 
weak interactions, such as \V c b\ and \ V u b\. As more data will be accumulated from the asymmetric 
.B-factories in near future, the theoretical and experimental studies on 5-meson decays would also 
give better understandings on the Standard Model and its possible extensions. 

In exclusive weak decay processes of hadrons, the effects of strong interaction are encoded in 
hadronic form factors. These decay form factors are Lorentz invariant functions which depend on 
the momentum transfer q 2 , and their behaviors with varying q 2 are dominated by non-perturbative 
effects of QCD. 

Over the past few years, a great progress has been achieved in our understanding of the exclusive 
semileptonic decays of heavy flavors to heavy flavors . In the limit where the mass of the heavy 
quark is taken to infinity, its strong interactions become independent of its mass and spin, and 
depend only on its velocity. This provides a new SU(2Nf) spin-flavor symmetry, which is not 
manifest in the theory of QCD. However, this new symmetry has been made explicit in a framework 
of the heavy quark effective theory (HQET) [Q]. In practice, the HQET and this new symmetry 
relate all the hadronic matrix elements of B — > D and B — > D* semileptonic decays, and all 
the form factors can be reduced to a single universal function, the so-called Isgur-Wise function 
[§, H|, which represents the common non-perturbative dynamics of weak decays of heavy mesons. 
However, the HQET cannot predict the values of the Isgur-Wise function over the whole q 2 range, 
though the normalization of the Isgur-Wise function is precisely known in the zero recoil limit. 
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Hence the extrapolation of q 2 dependences of the Isgur-Wise function and of all form factors is 
still model dependent and the source of uncertainties in any theoretical studies. Therefore, it is 
strongly recommended to determine hadronic form factors of B-meson decay more reliably, when 
we think of their importance in theoretical and experimental analyses. 

In this paper we develop a new approach for exclusive semileptonic B decays to D, D*, and 
predict the q 2 dependences of all form factors, inspired by the parton model. Previously the parton 
model approach has been established to describe inclusive semileptonic B decays 0, 3, and found 
to give excellent agreements with experiments for electron energy spectrum at all energies. The 
parton model approach for the inclusive decays pictures the mesonic decay as the decay of the 
partons in analogy to deep inelastic scattering process. The probability of finding a 6-quark in a B 
meson carrying a fraction x of the meson momentum in the infinite momentum frame is given by 
the distribution function f(x). 

While many attempts describing exclusive B decays often take the pole-dominance ansatze as 
behaviors of form factors with varying q 2 , in our approach they are derived by the kinematical 
relations between initial b quark and final c quark. According to the Wirbel et al. model ||, which 
is one of the most popular model to describe exclusive decays of B mesons, the hadronic form 
factors are related to the meson wavefunctions' overlap-integral in the infinite momentum frame, 
but in our model they are determined by integral of the distribution functions of quarks inside the 
mesons. 

In Section II, we develop the parton model approach for exclusive semileptonic decays of B 
meson. All the theoretical details for B — > Dlv and B — > D*lu are given in Section III. Section IV 
contains discussions and conclusions of this paper. 



2 Parton Model Approach for Exclusive Decays of B Meson 

We now develop the parton model approach for exclusive semileptonic decays of B meson by 
extending the previously established inclusive parton model, and by combining with the results of 
the HQET. Theoretical formulation of this approach is, in a sence, related to Drell-Yan process, 
while the parton model of inclusive B decays is motivated by deep inelastic scattering process. And 
the bound state effects of exclusive B decays are encoded into the hadronic distribution functions 
of partons inside an initial B meson and of partons of a final state resonance hadron. In Fig. 1, we 
show the schematic diagrams of Drell-Yan process and the related exlusive semileptonic decay of 
B meson. We assume that the Lorentz invariant hadronic decay width can be obtained using the 
structure functions, 

E B ■ dT(B -» D(D*)ev) = J dx J dy f B (x) E b ■ dT(b -» cev) f D (y) (1) 

following the Drell-Yan case. The first integral represents the effects of motion of b quark within B 
meson and the second integral those of c quark within D meson. The variables x and y are defined 
as fractions of momenta of partons to momenta of mesons, 

Pb = xp B , Pc = VPD , (2) 

in the infinite momentum frame, |p| 3> m. In this frame, the eq. (2) is valid for the four momenta 
of massive partons. The functions Jb(x) and /d(2/) are the distribution function of b quark inside B 
meson, and the fragmentation function of c quark to D meson, respectively, which are also defined 
in that frame. Since the momentum fractions and the distribution functions are all defined in the 
infinite momentum frame, we have to calculate the Lorentz invariant quantity, E ■ dT as defined 
in Eq. (1), to use at any other frame safely. 
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The distribution function can be identified with the fragmentation function for a fast moving 
6-quark to hadronize into a B meson in the infinite momentum frame. In general the distribution 
and fragmentation functions of a heavy quark (Q = t, b, c) in a heavy meson (Qq), which are closely 
related by a time reversal transformation, are of very similar functional forms, and peak both at 
large value of x. Brodsky et al. (jj have calculated the distribution function of a heavy quark, which 
has the same form as Peterson's fragmentation function ||. Therefore, here we follow the previous 
work of Paschos et al. Q to use the Peterson's fragmentation function for both distributions, 
and /r>(y). It has the functional form: 

f Q (z)=N Q z- 1 [l- 1 ---^-^ 2 , (3) 

where Nq is a normalization constant, and Q denotes b or c quark. This functional form is not 
purely ad hoc, but motivated by general theoretical arguments, that the transition amplitude for 
a fast moving heavy quark Q to fragment into a heavy meson (Qq) is proportional to the inverse 
of the energy transfer AE . 

This function (3) has a peak at z ~ 1 — y / eQ. Since the parameter e<g is related to the ratio of 
the effective light quark mass to the heavy quark mass ~ m 2 /mQ, the peak approaches to 1 in the 
limit of niQ — > oo, and the function shows the similar behavior to 5(1 — z). Hence we find that this 
functional form is supported by the HQET in this indirect manner. 

In the Drell-Yan process, the rest degrees of freedom of initial nucleons which do not take part 
in the scattering make incoherent final states, see in Fig. 1 (a). In the exclusive semileptonic 
decay of a heavy meson into a final state heavy meson, however, two sets of left-over light-degrees 
of freedom are summed to have the connection, i.e. see Fig. 2 (b). To connect them we need a 
relation between x and y from the decay kinematics. The momentum transfer of the decay between 
mesons is defined as 

q = PB~PD ■ (4) 
On the other hand, the momentum transfer of the partonic subprocess is given by 

= XPB - VPD ■ (5) 

Note that the second line of the Eq. (5) holds only in the infinite momentum frame while Eq. (4) 
holds at any frames. In fact, the heavy meson's momentum would be pn = Pq + k + 0(l/mg), 
where H = B,D or D*, and Q = b or c. And k denotes the momentum of the light-degrees of 
freedom, and has the size of the effective mass of a common light degree of freedom, A. Therefore 
we have q = g(p arton > up to the common part of the I/thq corrections in the infinite momentum 
frame, and with these kinematic relations, (4) and (5), we derive the following relation 
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y(x, q 2 ) = J x(x — l)m 2 B + (1 — x)q 2 + xm? D . (6) 

mc v 

Note that if x = 1, y = 1, this equation describes the decay of the free quark with mass to the 
final quark with mass mp. Substituting y of Eq. (1) for y(x,q 2 ) of Eq. (6), the double integral of 
Eq. (1) is reduced to the single integral over x. Using this relation, we can sum the intermediate 
states in Fig. 1(b). 

We note here that the kinematic relation (6) are valid approximations for the heavy-to-heavy 
resonance decays, with the common light-degrees of freedom of the size 0(\/rriQ). As explained 
before, in the limit where /q(x) = 5(1 — x) by increasing toq to infinity, we can reproduce the 
HQET leading term. By comparison, the inclusive parton model approach is more reliable for the 
heavy-to- light non-resonant decays to final states of many particles. 
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3 Form Factors for semileptonic B — > D(D*) Decays 
3.1 B -> Dez/ 

From Lorentz invariance we write the matrix element of the decay B — ► DeP in the form 

< J D|J^|i?>=/ + (g 2 )(p B +p D ) M + /_(g 2 )(pB-p D ^ , (7) 

and in terms of the HQET 

<D{v , )\J^\B(v)>=^^(i + (v-v')(v + v')^ + ^{v-v'){v-v')^ . (8) 

Due to the conservation of leptonic currents, the form factors multiplicated by do not contribute. 
Then the hadronic tensor is given by 

H„ v = < D\J^\B >< B\4\D > 

= ^\f+{<f)\ 2 {PB^PD v +PB v PD iJ ) , (9) 

and can be expressed by the Isgur-Wise function, 



(10) 



= R^lSiv ■ v')\ 2 {p B ^ Dv +p BvPD ^ ( 1 + 0(— ) ) , 

\ m Q / 

where 

_ 2^/m B m D 
mB + m D 

The partonic subprocess decay width for B — ► Dezv decay is given by 

E b ■ dF(b - cei/) = i(27r) 4 <5 4 (p 6 - p c - (z (parton) ) • 2G F 2 |y c6 | 2 ^ arton) ^ 

<i 3 p c d 3 p e d?p„ 



(2ir) 3 2E c (2n) 3 2E e (2ir) 3 2E u ' 
where 

V™ = 2 (p£j£ + Plpi ~ gTPe ■ Pu + ie^PeaPvp) , 

HftT* = N (p bfM Pc u +p bv Pc^) ■ (12) 

Here Bj^ Iton ^ denotes the partonic subprocess' hadronic tensor contributing B — > D decay. Now 
we need to discuss about it. In the inclusive decays B — > X c eu, the hadronic tensor is given by 

= 2 (i^Pc;, + PbvPcp ~ Q^uPb ■ Pc ~ ie^uafSPbPc) ■ 

As can be easily seen, this hadronic tensor contains all the possible spin configurations of b — > c 
transition from spin state. Here we are interested in only B — > D process, where the spin does not 
change during the process. Therefore, we have to choose only the spin-inert part which contributes 
to B — > D process from the inclusive hadronic tensor, W^ v . By comparing with Eq. (10), we find 
that the spin-inert part has the form of {pb^Pcy+PbuPc^)-, as shown in (12). Besides, we do not know 
how much spin-inert part really contributes to B — > D semileptonic decay, because other decays 
such as B — > D*, B — ► D** also contain spin- inert parts. Therefore, the parameter is introduced 
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to estimate the size of spin-inert part which contributes to B — > D process. The constant N will 
be later determined by the zero recoil limit of the Isgur-Wise function. 

Using the relation (2), we can write the hadronic tensor in the parton level as follows 

= Nxy (PB^PDv+PBvPDh) ■ (13) 

The momentum conservation of the partonic subprocess corresponds to the momentum conservation 
in the hadronic level in our model, as explained before. So we can substitute the Dirac delta function 
5 4 (pb — p c — gi(p arton )) for 5 4 (pb — Pd — o) m Eq. (11) with no loss of generality. Therefore, we write 
the decay width of B — > Dev, 

E B ■ dT(B -> Dev) = J dx f B (x) f D (y(x,q 2 ))E b dT(b -> cev) 

dx fB(x)fD(y(x,q 2 )) (2vr) 4 (5 4 (p B - p D - q) 

x G F 2 \V cb \ 2 N xy(x,q ) (pBuPDv + Pb u PDu)L^ u 
2/ 2x d 3 p D d 3 p e d 3 p v 
V [X,q ' {2Txf2E D (2iTf2E e (2TTf2E v ' 1 ' 

Hence we find the hadronic tensor 

HpviB -> Dev) = N J dx f B (x)f D (y(x, q 2 )) xy 3 (x,q 2 ) (p Btl PD u + Pb v Pd^ 

= N F(q 2 )(p B ^p Du +p Bu p D ^) , (15) 
where we factorize the function J-(q 2 ) defined as 

Hq 2 ) = [dx f B (x)f D (y(x,q 2 )) xy 3 (x,q 2 ) . (16) 



For given q 2 in our parton picture, the function J-(q 2 ) measures the weighted transition amplitude, 
which is explicitly related to the overlap integral of distribution functions of initial and final state 
hadrons. 

Comparing (15) with the Eq. (10), the Isgur-Wise function is calculated within the parton 
model approach 

\£(v-v')\ 2 (l + 0(— ) ) =4N-R-F(v-v') . (17) 
V m Q J 

As explained before, in our model the non-perturbative QCD effects are included through the 
distribution functions, so the function J-(q 2 ) contains all orders of 1/toq corrections. Note that, 
however, our model ignores the effects of the transeverse quark motion and off-mass shell, the full 
1/mQ effects are not being considered. From the value of the zero recoil limit of the Isgur-Wise 
function with 1/mq corrections ||, we can determine the numerical value of N. In Table 1, we 
show the numerical values of N with varying the parameters of the fragmentation functions, (c^, e c ). 
Unfortunately we cannot obtain the analytic form of J'iq 2 ) due to the nontriviality of the integrals 
of f(x). However, we can numerically obtain the behaviors of the Isgur-Wise function with varying 
q 2 from the Eq. (17), and calculate its slope parameter. The normalization constant N and the 
slope parameters p 2 with the input values of (e&, e c ) are also shown in the Table 1. We varied e\, 
from 0.004 to 0.006, and e c from 0.04 to 0.1. As explained earlier, the parameter eg is related to 
the ratio of the effective light quark mass to the heavy quark mass ~ m 2 /mQ. And therefore, our 
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input values of (e&, e c ) correspond to the range m 2 /m 2 = 0.04 ~ 0.15, which is consistent with the 
prediction of the HQET Q. For more details on the comparison with the HQET, see Section IV. 
The q 2 spectrum is given by 

dT(B^Deu) _ G F 2 \V cb \ 2 2 / 2 2 2 2 2 2 n3/2 



967r 3 m^ 



^(g 2 )iV ((m| - m 2 , + g 2 ) 2 - 4m|g 2 ) . (18) 



And in Fig. 2(a), our prediction is shown with the parameters (e& =0.004, e c =0.04), compared 
with Wirbel et aVs model prediction ||, which shows quite a good agreement. 

3.2 B -> D*eu 

We write the matrix element of the decay B — > D*ei> in familiar form 

< D*\V, + A,\B> = - * 

ma + rriD* 

+(m B +m D *)e* ^Ai(q 2 ) ■ (p B + PD*)^A 2 (q 2 ) 

-2m D *^^A 3 (9 2 ) + 2m D .^g M A (g 2 ) , (19) 



where 



and in terms of the HQET 

<D*(i/)\V lt + A ll \B(v)> = y/m B m D . 

X(itiv( v ■ v')e liVafi e fV v' a v^v{q 2 ) + £ M (v • "Ofr " ^ + 

-U^-w'K-w^-UCf^e*-!;^) . (20) 

In fact, A()(q 2 ) and A^(q 2 ) do not contribute here because of leptonic currents conservation. With 
the heavy quark expansion, the hadronic form factors are related to the Isgur-Wise function £(v-v'), 
such as £i(v ■ v') = £(v ■ v')(ai + C(l/mg)), where «a 2 = and aii = 1 otherwise. 
The hadronic tensor of B — > D*ev decay is obtained using the HQET 

R »» = J2 <d *\ v » + a k\ b><d *\ v » + a "\ b> * ( 21 ) 



R*- l \t(vv')\ 2 



{m B + m D *y p p Mq 



where i?* = 2^m B m B * /{m B + m B *). Investigating the above relation, we can find the form of 
the hadronic tensor for the partonic subprocess which contributes to the decay B — ► D*eu. The 
resulting hadronic tensor is written down in the form, 



/ 2a 2 

^(partem) / ^ 

I (m/i I 111 it- )' 2 



N l {Pb^Pcy + PbvPc 

m B + m D *) 2 



2 \ 1 " 7 r \2 1 (^2 Wfe • - iN 3 e^afsp^pP) . (22) 
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The parameters iVj's give the relative size of form factors and overall normalization. In general they 
are not constants and have the q 2 dependences. In the heavy quark limit, iVj's become constants 
and the values are equal to that of the normalization constant N in B — > Dev process. 

In order to investigate the procedure more conveniently, we define the ratios of form factors as 
follows: 

I ^ 



Ri = 



R2 



V(q 2 



(m B + m D * 



1 



Ai(q 2 ) 

q 2 \ A 2 (q 2 ) 

(m B + m D *) 2 ) Ai(q 2 ) 



(23) 



where V{q 2 ), A\(q 2 ) and A 2 (q 2 ) denote vector and axial vector form factors respectively. Then we 
can write the relations among form factors and the Isgur-Wise function as 

- (l 



Mr) 

A 2 (q 2 ) 
V(q 2 ) 



(m B + m D * 
R2R*~\{q 2 ) , 

RiR*-^(q 2 ) ■ 



R*-^(q 2 



(24) 



Note that the Isgur-Wise function ^(v ■ v') in these expressions contains full 1/rriQ corrections and 
it corresponds to £(u • v') in the Ref. ||, which is normalized at zero recoil up to the second order 
corrections £(1) = 1 + 8i/ mQ . And R± and R 2 become to be unity in the heavy quark limit, and 
Neubert's estimates of the 1/rriQ corrections Q give 

i?i » 1.3 , R 2 ~ 0.8 , (25) 



which are model dependent. Recently CLEO |l0| obtained the values of the parameters R%, R 2 by 
fitting them with the slope parameter of the Isgur-Wise function p 2 . Since Ri and R 2 have very 
weak q 2 dependence, in the CLEO analysis they approximately obtained the q 2 independent values, 



Ri = 1.30 ±0.36 ±0.16 
R 2 = 0.64 ±0.26 ±0.12 



(26) 



Hereafter we also take them to be constants for simplicity. 

In our model the parameters JVi, and ./V3 are represented in terms of i?i and R 2 as follows, 

„2 



N/2 



+ 



N 7 



1 — 2q 2 1 (tub + mD*) 2 

\ £ )( 

(m B + m D *) 2 ) \ 
A + R\) + 



(m B + m D *) 2 



■2R\ 



(1+R 2 ) 2 m 2 B + q 



N 
~2 
NRi 



2m B m B) * 



m B + m 



+ 



r(l 



2m 2 D , 



-{l-R 2 f + 



2m B mD* — q 
2m 2 D , 



(1 



R 2 ) 



Ri) 



D* 



In the heavy quark limit, we know that R\ = R 2 = 1. Using the expression Ri 
can separate the leading contributions and 1/mQ corrections in AVs: 



(27) 

1 ± 0(l/m Q ) we 



Nl 

N 2 
N 3 



N 

T (l±^ 2 ) 



+ 0(1 



E(l + R 2 ) + 0(1 
NRi . 



-R2) 
Ri) 



(28) 
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When Ri, R 2 -> 1, we explicitly see that Nx = N 2 = N 3 ^ N. 

Now substituting £(v ■ v') for J-(q 2 ) with the Eq. (17), the hadronic tensor for the decay 
B —* D*ev is given by 



^ 2 ) 
2(1- 



(1 



2<? 2 



(ms + m D .) 2 



)Nl(pBuPD* v +PB v PD* 



(m B + m D *) 2 



(29) 



The q 2 dependences of form factors are mainly determined by the function ^(q 2 ), instead of the 
commonly used pole-dominance ansatze. 

When we calculate the hadronic tensor within the HQET framework, we have generally some 
parameters parametrizing non-perturbative effects, which are obtained in model dependent ways. 
The slope parameter is such a characteristic parameter of the Isgur-Wise function, which represents 
the common behaviors of form factors. We calculated it, and find that the value of the slope 
parameter is related to the parameters and e c in our approach. The HQET also contains the 
parameter Ai ~ ~(^q) which is related to the kinetic energy of the heavy quark inside the heavy 
meson pl[| , and spin-symmetry breaking term A2 = ji^n 2 / — m 2 p) corresponding to the mass splitting 
of pseudoscalar mesons and vector mesons. In our approach, we have two parameters R\ and 
R 2 correspondingly, which are experimentally measurable. Also the values of and e c can be 
independently determined from the various methods experimentally and theoretically. We use the 
values for R\ and R 2 from the CLEO fit results of Ref . flCfl , and for ej, and e c from Ref . [||, ^] . 

Finally, we obtain the decay spectrum 



dT(B -> D*eu) 
dq 2 



G 2 F \V r h \ 2 
192vr 3 



m 



Hq 2 



B 



2 1 2\2 

m D * +q ) 



m 2 B W±(q 2 ) ((m 

3 
2 



m 2 D * 



2\2 



+ q 2 ) 



, 2 2\ 1 ' 2 



2 2^ 

m B q ) 



A m2 B W ^1 2 ) « ~ ™l* + Q 2 ) + 3™ 2 B W 3 (q 2 ) 



(30) 



where 

W^q 2 ) 

W 2 (q 2 ) 

W 3 (q 2 ) 



-iVi 1 



2q 2 



(m B + m D *) 2 ^ 
Nl (m 2 B -m 2 D , +q 2 ) ( 1 - 



2q 2 



(m B + m D *) 2 



2N 3 q 2 1 



(m B + m D *) 2 



-N im 2 B q 2 I 1 



2q 2 



v (m B + m D *) 2 J 
+ N 3 q 2 {m 2 B -m 2 D , + q 2 ) 1 1 



{m B + m D *) 2 



+ N 2 q 2 (m 2 B +m z D ,-q z ) I 1 



2q 2 



(m B + m D *) 2 



(31) 



and J~(q 2 ) is defined in (16). The result is plotted in Fig. 2 (b), also compared with the CLEO 
data ]l0|| . The thick solid line is our model prediction with the parameters (e;, =0.004, e c =0.04) 
and the dashed line the Wirbel et al. model prediction B. 
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4 Discussions and Conclusions 



All form factors show the common behavior for varying q 2 , which is described by the Isgur-Wise 
function of the HQET, which represents the common non-perturbative dynamics of weak decays 



of heavy mesons. Ever since Fakirov and Stech 13], the nearest pole dominance has been usually 



adopted as the dependence of common behaviors on q . In our approach, their q dependences 
are extracted from the kinematic relations of b- and c-quark. When 6-quark decays to c-quark, 
the momentum transfer to leptonic sector is equal to the difference between 6-quark momentum 
and c-quark momentum in the parton picture. The b- and c-quark momenta within the B and 
D mesons have some specific distributions. For given momentum transfer q 2 , there exist possible 
configurations of b- and c-quark momentum pairs (pb, p c ), and each pair is appropriately weighted 
with the momentum distributions of the quarks. Our J~(q 2 ) function in (16) measures the weighted 
transition amplitude for given q 2 in the parton picture; it is explicitly given by the overlap integral 
of distribution functions of initial and final state hadrons. This is common to all form factors, as 
explained in Section II. 

As mentioned earlier, the non-perturbative strong interaction effects are also considered through 
the distribution functions in our model, so £(v ■ v') obtained from Eq. (17) corresponds to the 
hadronic form factor £(f • v') defined in the Ref. ||, including 1/nib corrections rather than the 
lowest order Isgur-Wise function. And the slope parameter of our results in Table 1 also corresponds 
to p 2 related to £(v ■ v'). We obtain the values of the slope parameter p within the parton model 
framework, as in Table 1, 

p 2 = 0.552 - 0.858 , 

which are compatible with the Neubert's prediction ||], 

p 2 ~ p 2 ±0.2 = 0.7 ±0.2 . 

Our result is rather smaller than the predictions of other models, 

p 2 = 1.29 ±0.28 [14] , 
p 2 = 0.99 ±0.04 [T 



but it is consistent with the average value measured by experiments [16] 



p 2 = 0.87 ±0.12 . 

In calculating the numerical values, we still have two free parameters e& and e c of Eq. (6), i.e. 
of the heavy quark fragmentation functions. Their values can be determined independently from 
the various experimental and theoretical methods^. For the parton model to be consistent with 
the HQET, we require that with the fixed value of the parameter e<g, all the appropriate results of 
the parton model approach agree with those of the HQET. In other words, the value of parameter 
€q should be determined to give all the phenomenological results to coincide with those of the 
HQET. In this point of view, we have previously studied the parton model approach for inclusive 
semileptonic decays of B meson in the Ref. |L7|], and showed that the value €b ~ 0.004 gives 
consistent results with those of the HQET. In this paper we use the value of as 0.004 or 0.006. 
The latter value is given by the experiments for the determination of the Peterson fragmentation 
function [12]. We find that our prediction of the slope parameter p 2 with the parameter e;, = 0.004 



5 Theoretically, it is interesting to calculate /q(z) = 8(1 — z) + 0(l/m 2 Q) within the HQET, which can reproduce 
the phenomenological predictions. 
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and e c = 0.04 gives the best agreed value ~ 0.7 with that of the HQET, p 2 = 0.7 ± 0.2. In this 
context, we conclude that our model with the parameter e& = 0.004 gives consistent predictions 
with the HQET. 

To investigate the possible model dependences of our approach on the input fragmentation func- 
tions, we now study the cases with other fragmentation functions. Following the Artru-Mennessier 
string model fli~8|| for heavy flavors, we have the fragmentation function such as 

(l-z) a ( bm 2 n \ 
/ Q (,)=iV^^ex P (--^j , (32) 

with parameters b = 0.8 GeV~ 2 and a ~ 0.5, and the c- and 6-quark masses identified with the 
masses of the lowest-lying vector meson states fL9| . Another fragmentation function is recently 
developed with the help of the perturbative QCD (PQCD). Inspired by the PQCD of the heavy 
quarkonium, the functional form for the heavy-light mesons |2(J is derived as 

r Q z(l-zf f co j,j 



Iq( z ) = N Q Q _ i _ -j-y 6 (6 - 18(1 - 2r Q )z + (21 - 74r Q + 68r 2 



z 



-2(l-r Q )(6-19r Q + 18r 2 Q )z 3 + 3(l-r Q ) 2 (l-2r Q + 2r 2 Q )z 4 ) , (33) 

with the free parameter tq being the mass relation mQ/(rriQ +m q ), and with the light quark mass 
m q . We call this function the PQCD inspired fragmentation function. With these two fragmentation 
functions, we obtain the numerical values of p 2 = 1.124 from the string fragmentation function, 
and p 2 = 0.605 from the PQCD inspired fragmentation function. The results are shown in Table 2 
compared with those with Peterson's function. As can be seen, the string fragmentation function 
gives a somewhat larger value of p, but the 'recently developed' PQCD inspired function gives a 
similar value to that of Peterson fragmentation function. 

Phenomenologically, our model prediction on q 2 spectrum in the B — > D*ev decay shows a good 
agreement with the result of the CLEO |l(J, as shown in Fig. 2(b). If we let /q(z) = 5(1 — z) and 
Rl = i?2 = 1 in our model, we can reproduce the lowest order results of the HQET, and obtain 
the similar plot with those of other models in Fig. 2(b). For the B — > Dev decay, our results agree 
with those of other models, as in Fig. 2(a). 

Finally we obtain the ratio of integrated total widths T(B — > D*)/T(B —* D) ~ 2.66, which 
agrees with the experimental results plf . It may be a phenomenological support of our model 
because this quantity is independent of the CKM elements |V^,|, which has uncertainties in de- 
termining its value yet. The perturbative QCD corrections can be factorized in the decay width 



calculation |], |22|, which does not affect the ratio T(B —* D*)/T(B — > D). 

To summarize, we investigated the form factors for exclusive semileptonic decays of S-meson 
to D, D*, based on the parton picture and helped by the results of the HQET. We obtained the 
numerical results for the slope of the Isgur-Wise function, which is consistent with the experimental 
results, and we extracted the dependences of hadronic form factors on q 2 by varying input heavy 
quark fragmentation function without the nearest pole dominance ansatze. 
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Table 1: The values of the normalization constant N and the slope parameter p 2 are shown with 
several choices for the input parameters e\, and e c . 







£fe = 


0.004 






£b = 


0.006 






e c = 0.04 


0.06 


0.08 


0.1 


0.04 


0.06 


0.08 


0.1 


N 


0.938 


1.068 


1.190 


1.306 


1.105 


1.223 


1.338 


1.449 


P 2 


0.705 


0.646 


0.609 


0.582 


0.896 


0.844 


0.810 


0.785 



Table 2: Comparison of the slope parameter p 2 for the different fragmentation functions 



Fragmentation Function (F. F.) slope parameter p 2 

Peterson F. F. 0.5521 - 0.8582 

String F. F. 1.1240 

PQCD inspired F. F. 0.6053 
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Figure 1: (a) The diagram of Drell-Yan Process, (b) The schematic diagram of semileptonic 
exclusive decay of 5-meson in the parton model. 



Figure 2: (a) q 2 spectrum in the B — > Dev decays. The solid line is our model prediction and the 
dotted line the Wirbel et al. model prediction ||. We used the values of parameters, (e& =0.004, 
e c =0.04). (b) q 2 spectrum in the B — > D*ev decays. The solid line is our model prediction and the 
dotted line the Wirbel et al. model prediction |J. We used the values of parameters, (e& =0.004, 
e c =0.04). The data are quoted from Ref. [ 10 1 . 
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Fig. 1 (a) 
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